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Abstract 

Extending the super symmetric method proposed by Tkachuk to the complex do- 
main, we obtain general expressions for superpotentials allowing generation of quasi- 
exactly solvable PT-symmetric potentials with two known real eigenvalues (the ground 
state and first-excited state energies). We construct examples, namely those of com- 
plexified non-polynomial oscillators and of a complexified hyperbolic potential, to 
demonstrate how our scheme works in practice. For the former we provide a connec- 
tion with the sl(2) method, illustrating the comparative advantages of the supersym- 
metric one. 



Running head: SUSY and PT-symmetric potentials 
PACS: 03.65.Fd, 03.65. Ge 

Keywords: supersymmetric quantum mechanics, quasi-exactly solvable potentials, PT sym- 
metry, non-polynomial oscillator, hyperbolic potential 

*Directeur de recherches FNRS 

1 



1 Introduction 



Non-Hermitian Hamiltonians, in particular the PT-symmetric ones, are of great current 
interest (see e.g. [l j jL^j and references quoted therein). The main reason for this is that 
PT invariance, in a number of cases, leads to energy eigenvalues that are real. In this 
regard, the early work of Bender and Boettcher Jl| is noteworthy since it has sparked off 
some very interesting developments subsequently. 

Some years ago, Tkachuk [0 proposed a supersymmetric (SUSY) method for generating 



quasi-exactly solvable (QES) potentials with two known eigenstates, which correspond to 
the wavefunctions of the ground state and the first excited state. A distinctive feature of 
this method is that in contrast with other ones, it does not require the knowledge of an 
initial QES potential for constructing a new one. Later on, the procedure was extended to 
deal with QES potentials with two arbitrary eigenstates |TJ| or with three eigenstates [ITH] . 



Quite recently, Brihaye et al. |16j established a connection between the Tkachuk approach 



and the Turbiner one, based upon the finite-dimensional representations of sl(2) |17| 



In this letter, we pursue Tkachuk's ideas further by considering an extension of his 
results to the case of PT-symmetric potentials. As a consequence, we arrive at a pair of 
solutions, one of which is obtained by a straightforward and natural complexification of 
Tkachuk's results while the other is new. We demonstrate the applicability of our scheme 
by focussing on two specific potentials, both of which are PT-symmetric. 

2 Procedure 

2.1 The Hermitian case 

Let us start with the Hermitian SUSY case, where the two known eigenstates are the ground 



state and the first excited state. As is well known ||18|| , the SUSY partner Hamiltonians are 
given by 

#(+) = AA = — f- + V {+ \x), ifH = AA= — f- + v(-\x), (1) 
ax 2 ax z 
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for a vanishing factorization energy. Here A and A are taken to be first-derivative differential 
operators, namely 

A=^- + W(x) } A=—^ + W(x), (2) 
ax ax 

where W(x) is the underlying superpotential and V^- ± \x) are the usual SUSY partner 
potentials 

V {±) (x) = W 2 (x) =F W\x). (3) 

We assume SUSY to be unbroken with the ground state of the SUSY Hamiltonian 
H s = diag(#W,#H) to belong to ffW; 

H<- + Vo +) (x) = 0, 4 +) (x) = C { +) exp (- J X W(t)dt) , (4) 

Cq being the normalization constant. Equation @ implies 

sgn(TU(±oo)) = ±1. (5) 

In the above formulation of SUSY, the eigenvalues of and are related as 

E^^E^, ^S +) =0, (6) 

while the corresponding eigenf unctions are intertwined according to 

Following Tkachuk let us consider expressing if* - ) in the form 

#H = + e = A x Ai + e, (8) 

where e = i?i = Eq corresponds to the energy of the first excited state of or of the 
ground state of H^' and the operators A\ and A\ are such that relations analogous to (0) 
hold in terms of a new superpotential W\(x). Thus we can write 

V { -\x) = v} + \x)+e (9) 

with = W?(x) - Wftz) from (H). 
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The ground state wave function of H[ + ^ (or )) can be read off from 

= Ct ] exp (- W!(t)dt 



as 



(10) 



where Cq "* is the normalization constant and 



sgn(Wi(±oo)) = ±1. 



ti: 



On the other hand, the wave function of the first excited state of is obtained in the 
form ^ ^ 

4 + \x) = Cf^exp (- f W 1 (t)d?j , c{ +) = ^L-, (12) 

as guided by (|7|) and fllPf) . 

Using (|3|) and (||), it is clear that the superpotentials W(x) and W\{x) have to satisfy 
a constraint 

TU 2 (x) + W(z) = Wl{x) - W[{x) + e. (13) 

The problem therefore amounts to finding a set of functions W(x) and Wi(x) satisfying 
fll^D for some e > 0, along with the conditions (|5]) and ([11]). For this purpose, Tkachuk 
introduced the following combinations 



W±(x) = W 1 (x)±W(x), 

which transform ([H|) into 

W' + {x) = W-(x)W+(x) + e. 
An advantage with Eq. (|T3j) is that it readily gives in terms of W + {x): 



Wjx) 



Consequently the representations 



W(x) 



x 



W' + {x) - £ 

W + (x) 



W' + {x)-e 
WJx) 



WAx) 



WJx) + 



W' + (x)-e 
W+(x) 



(14) 



(15) 



(16) 



(17) 



satisfy Eq. (|i~3|). In fllTD , W+(a;) is some function for which VF(x) and Wi(x) fulfil condi- 
tions (|5|) and (|TT]) . This means 



sgn(W+(±oo)) = ±1. 



Restricting to continuous functions W+(x), the above condition reflects that W+(x) 
must have at least one zero. Then from ( ]T5| ) and (|T7|), W-(x), W(x), and Wi(x) may have 
poles. Tkachuk considered the case when W + (x) has only one simple zero at x — Xq. In 
the neighbourhood of xq, one gets 

W' + {x) - e W' + {x ) -e + {x- x )W^(x ) + ■■■ 



W+(x) (x - X )W{{x ) + • • • 

so that the superpotentials will be free of singularity if one chooses 



(19) 



e = W' + {x ). (20) 

One is then led to 

W' (>A - W (<r„\~\ 

WAx) 



W(x) 9 



1 w , x) Wj(x)-Wj(x o y 1 



WAx) 



Wjx) 



2 

"(21) 

To summarize, provided the continuous function W+(x) with a single pole at x = x is 
such that W(x) and Wi(x), given by (pip, satisfy conditions (0) and (O), the Hamiltonian 
has two known eigenvalues and e, given by fl20D , with corresponding eigenfunc- 
tions @ and 

ij[ + \x) = c[ +) W + (x) exp (- J X Wi(t)dtj . (22) 
In deriving fl2"2|), use is made of (0), (|T2"D, and (0). 

2.2 The non- Hermit ian case 

In the non-Hermitian case, we have to deal with complex potentials. Consider the decom- 
position J| 

W(x) = f(x) + ig(x), V i±] (x) = V^\x) + iV} (±) (x), (23) 
where /, g, Vj (±) e R and 

v£ ] = f-g 2 Tf\ v} ±) = 2fg T g'. (24) 
If V^'(x) are PT-symmetric, then f(x) and g(x) are odd and even functions, respectively. 
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It may be noted that Eqs. (^) - (|T3| ) remain true for some real and positive e. Employing 
the separation 



W 1 (x) = f 1 (x)+ig 1 (x), v} +) (x) = V}p (x) + iV}+ ] (x) 



(25) 



with fx, gi, V X j G R, the behaviour of fi(x) and gi(x) also turns out to be odd and 

even, respectively, should V' + ^(x) be PT-symmetric. 

In the non-Hermitian case, the conditions (f|) and are to be replaced by 



sgn(/(±oo)) = Bgn(/i(±oo)) = ±1. 



(26) 



These conditions are actually compatible with the odd character of / and fx. 

On introducing the first relations of ( p3|) and ( p5[ ) into (0) and splitting into real and 
imaginary parts, we get the system of two equations 



f -9 2 + f' = fi-9 2 i-f[ + e, 2fg + g' = 2f 1 g 1 -g' 1 . 



(27) 



The problem now amounts to finding a pair of odd functions /, fi and a pair of even 
functions g, gi, satisfying Eq. fl2T|) for some e > 0, as well as the conditions (|26|). 
To this end, we introduce the linear combinations 



f±{x) = fi{x) ± f{x), g ± (x) = gi{x) ± g(x), 
which replace the constraints (|27| ) by 

/+ = /_/+ - g-g+ + e, g' + = f+g- + /_#+. 
Solving for /_ and g_ in terms of /+ and g + , we get 



(2f 



(29) 



/- 

Hence the functions 

fi = \ /+ + 



(/;- 


- e)f+ + g'+g+ 




fl + 9l ' 


(/; 


- e)f+ + g'+g+ 




fl + 9l 


(/; 


- e)f + + g' + g + 



if l + -e)g + + g' + U 

ft + gl 



(30) 



n+9i 



9= 2 



9i 



9+ ~ 
9+ 



(/; - 1)9 + - g' + h 

fl + 9l 

(/; - ^+ - 4/j 



(31) 



satisfy the coupled equations (|27|) . In (|3T|), /+ must be such that the conditions (^) are 
fulfilled. These suggest 

sgn(/ + (±oo)) = ±1. (32) 

If we restrict ourselves to continuous functions f+(x) and g+(x), the condition (|32]) 
shows that f+(x) must have at least one zero. For simplicity's sake, we assume that f+(x) 
has only one simple zero at x = Xq. This means that the parity operation is defined with 
respect to a mirror placed at x = xq. Thus in the neighbourhood of xq, we get 

(f + -e)f + + g'g + 



t \ [f+M - e]/+(x ) + g + (x )g'i(x ) 

[X — X ) ; — h ■ 



(33) 



if g+{x ) ^ 



+ 



£ — f+[ x o) 
U+ - e]g+ - g' + U 

fl + gl 

_ /;(x )-e | 

g+( x o) 

[f + (x ) + e]gUx ) 



if fif+(x ) = 0, 



(34) 
(35) 



if g+(x ) ^ 
H if #+(x ) = 0. 



(36) 
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Hence the superpotentials will be free of singularity if either g+{ x o) ^ and e arbitrary or 
g+{x ) = and e = /+(%) (or e = H^Zo) since S-+(x ) = 0). 
The superpotentials may therefore be written as 



W(x) 



W+(x) - 



W'Ax) 



W + {x) 

if W + (x ) = ' i g+( x o) 7^ and e > 0, or 



WAx) 



W+(x) + 



W'Ax) 



WAx) 



(37) 



W(x) 



W+(x) 



W + {x)-W' + {x Q ) 
W + {x) 



Wxix) 



W+(x) + 



W[(x)-W[(x ) 
WAx) 



(3? 



if W+(x ) = ig+{x ) = and e = W' + (xq). 

To summarize, in the non-Hermitian PT-symmetric case, we get two types of solutions: 
the one given by fl38|) is obtained as a straightforward consequence of the complexification 
of Tkachuk's result, while the other given by fl3~7|) is new. 



3 Applications 

3.1 A family of complexified non-polynomial oscillators 

As the first application of our scheme we consider the case of a family of complexified non- 
polynomial oscillators. This corresponds to the choice /+ = ax, g + = bx 2m (a > 0, b ^ 0, 
m EN), which leads to 



W+(x) = ax + ibx 



2m 



m G N. 



(39) 



Here x Q = and the first type of solutions, given in ([$7]), applies to the case m — 0, while 
the second one, given in (|38|) , has to be used for m G No. 

The former case reduces to the well-studied exactly solvable PT-symmetric oscillator 
potential ||, [II]. O], thus partly accounting for the name of the family of potentials. Indeed, 
to get it in the standard form we have to set a = 2, b = —2c, and e = 4a. Then 
assumes the form 



V i+ \x) = (x-ic) 2 + 2(a-l) + 



(x — icY 



along with 



W(x) 



x 



ic + 



a 



WAx) 



a — 



x — 1C 



(40) 



(41) 



X — 1C x — 1C 

These agree with the two independent forms of the complex superpotentials proposed by 
us previously |l l]] in connection with para-SUSY and second-derivative SUSY of (|40|) . The 
potential (|40D can be looked upon as a transformed three-dimensional radial oscillator for 
the complex shift x — > x — ic, c > 0, and replacing the angular momentum parameter / by 
a — h. The presence of a centrifugal-like core notwithstanding, the shift of the singularity 
off the integration path makes flip]) exactly solvable on the entire real line for any a > 
like the harmonic oscillator to which ( |40| ) reduces for a = | and c = 0. 

In contrast, the case m G No is entirely new. For such m values, the superpotentials are 
given by 



W(x) = ~ 



ax + ibx 2m - 



2mibx 



2m-2 



a + ibx 2m 1 
where we used e = WL(0) = a. 



W x {x) = - 



ax + ibx 2m + 



2mibx 



2m-2 



a + ibx 



2m-l 



(42) 



The corresponding potentials turn out to be 



x 



1 

4 

+ 



-b 2 x im + 2iabx 2m+1 - 8mibx 2m ~ l + a 2 x 2 - 2a + 



Am{m — l)ibx 2m 3 



a + ibx 



2m-l 



4m(m — l)iabx 



2m-3 



m = 1,2, 3, 



(43) 



(a + ibx 2 ™- 1 ) 2 

which are QES and are seen to be PT-symmetric as well. The ground and first-excited 
state wave functions corresponding to the above QES potentials are easily determined to 
be 



.1 2 _ ^ x *rn + l 

4 2(2m + l) 



oc (a + ibx 2 ™- 1 )™^ 2 " 1 ^ exp 

i/j{ + \x) oc x{a + ibx 2 " 1 - 1 )^ 1 - 1 ^^ 1 -^ exp 

It is worth noting that the first member of the set (^) obtained for m = 1 

i 

V {+ \x) = - (-&V + 2iabx 3 + a 2 x 2 - 8ibx - 2a 



4 2(2m+l) 



(44) 
(45) 

(46) 

is a quartic potential differing from the known QES ones 0, || . All the remaining members 
of the set, starting with that associated with m = 2, 



.r 



-&V + 2iabx 5 



I6ibx 3 + a 2 x 2 -2a + 



8ibx 



+ 



8iabx 



(47) 



a + ibx 3 (a + ibx 3 ) 2 

are non-polynomial potentials. As for the PT-symmetric oscillator (f40D , the shift of the 
singularity off the integration path makes such potentials QES. 

On introducing the new variable z = x(a + ibx 2m ~ 1 )~ 1 ^ 2m ~ 1 \ the first-excited state wave 
function ( f4"ol) can be rewritten in terms of the ground state one fl4*4] ) as ip[ + \z) oc zip^ r \z). 
By setting in general ip^(z) = ip^ (z)^W(z), where 0o + ^(,2) oc 1 and (f>i\z) oc 2, the 
Schrodinger equation for the potentials (|43|) is transformed into the differential equation 



T<f>M{z) = 



_ a -2/(2m-l) (1 _ ihz 2m-lYm/{2m-l)^ + ^ 



0i +) W = ^ +) 0i +) W- (48) 



For m = 1, the coefficient of the second-order differential operator in ([48D becomes a 
quartic polynomial in z, thus showing that T can be expressed as a quadratic combination 
of the three sl(2) generators 

o d , d N d 



dz dz 2 dz 



(49) 



corresponding to the two-dimensional irreducible representation (i.e., iV = 1 in fl49|)) fT6 
P!7[ . The result reads 

T = cT 2 (~b 4 Jl - 4ifo 3 J+ J + 6b 2 J + J_ + 4i6J J- - J 2 - 2i& 3 J+ + 6b 2 J + 2ibJ_ + 36 1 



(50) 

For higher m values, the differential operator T contains a non-vanishing element of 



the kernel [16|. It is worth stressing that in such a case, the sl(2) method becomes quite 
ineffective for constructing new QES potentials, whereas the SUSY one is not subject to 
such restrictions. 

3.2 A complexified hyperbolic potential 

Our next example is that of a complexified hyperbolic potential induced by the represen- 
tations f + = A sinh ax, g + = B (A, a > 0, B ^ 0). We then get 



W+(x) = Asinhax + iB, 



(51) 



which gives for x = 

W(x 
WAx 



A sinh ax + iB — 
A sinh ax + iB + 



Aa cosh ax — e 
A sinh ax + iB 
Aa cosh ax — e 



A sinh ax + iB 

The resulting expression for the complexified hyperbolic potential is 

fx) = - A 2 sinh 2 ax - 4Aa cosh ax + 2e + a 2 - B 2 + 2iAB sinh ax 

e 2 _ a 2 (j4 2 _ £2)- 



(52) 
(53) 



+ 



(54) 



(A sinh ax + i£>) 2 

Clearly V™ (x) is PT-symmetric. The two known eigenstates of fl51|) correspond to the 
ground state and first excited state as outlined earlier. These are 

ifjQ + \x) oc (Acoshax — z/)*\ "^/(Acoshax + v)ii l+ ^) 

'A 1 i M . , 

x exp cosh ax iBx arctan — smh ax 

F 2a 2 2 V5 
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oc 



le / v 

H arctan — tanh ax 

2av \B 

ax {A 

oc cosh — exp cosh ax 

2 P V la 



if < B 2 < A 2 



if B = 0, 



V cosh ax exp 
if £? 2 = A 2 , 
oc (_B 2 + A 2 sinh 2 ax) 1/4 exp 



^4 1 e 

cosh ax i(5v4x H (sech ax + id tanh ax) 

2a 2 2Aa v 7 



cosh ax - 

2a 2a/i 



e / A cosh ax 
arctan 



a 



1 i /A , \ ie , fa 

iBx arctan — smh ax H arctanh — tanh ax 

2 2 \B J 2a/i 

if B 2 > A 2 , 



(55) 



and 



ijj[ + \x) oc ( A sinh ax + iB) (A cosh ax — z/) 4 ( x <*" ) ( A cosh ax + v) 
x exp 



A l._ i M . , 

cosh ax ii)X H — arctan — smh ax 

2a 2 2 \B 



ie / v 

H arctan — tanh ax 

2av \B 

, ax (A 

oc smh — exp cosh ax 

2 P V 2a 



if < B 2 < A 2 , 
if B = 0, 



oc (sinh ax + i<5)"\/sech ax 
x exp 



1 e 

cosh ax iSAx ; — (sech ax + id tanh ax) 

2a 2 2Aa v y 



if B 2 = A 2 



oc (A sinh ax + LB) (5 2 + A 2 sinh 2 ax) 1/4 



x exp 



A e /y4coshax\ 1 

cosh ax H arctan ii5x 

2a 2a// V /x / 2 



H — arctan — smh ax 
2 \B 



ie / /i 

arctanh — tanh ax 

2a/x 



if B 2 > A 2 , (56) 



where i/ = \/ A 2 — B 2 , /j = \J B 2 — A 2 , and 5 = sgn(S). 

In (ES|) and fl5D|), we have included the case B = for which the potential V^ + '(x) of 
Eq. (|54D reduces to one of the potentials studied by Tkachuk fl3|| , which itself is a special 
case of the Razavy potential |19| . 
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4 Conclusion 



To conclude, we have carried out in this paper a complexification of the SUSY method 
proposed recently by Tkachuk. This allows us to generate QES PT-symmetric potentials 
with two known real eigenvalues. We have also constructed suitable examples, namely those 
of a family of complexified non-polynomial oscillators and of a complexified hyperbolic 
potential, which serve to illustrate the viability of our scheme. For the former, we have also 
provided a connection with the sl(2) method, which illustrates the comparative advantages 
of the SUSY one. 
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